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We examine numerically the process ol gravitational wave recoil in the merger of two black holes 
in non head-on collision, in the realm of Robinson- Trautman spacetimes. Characteristic initial data 
for the system are constructed, and the presence of a common apparent horizon implies that the 
evolution covers the post-merger phase up to the final configuration of the remnant black hole. The 
net momentum flux carried out by gravitational waves and the associated impulses are evaluated. 
Our analysis is based on the Bondi-Sachs conservation laws for the energy momentum of the system. 
The net kick velocity Vk imparted to the merged system by the total gravitational wave impulse is 
also evaluated. Typically for a non head-on collision the net momentum flux carried out by gravi- 
tational waves is nonzero for equal-mass colliding black holes. The distribution of Vk as a function 
of the symmetric mass ratio n is well fitted by a modified Fitchett ^-scaling law, the additional 
parameter modifying the law being a measure of the nonzero gravitational wave momentum flux 
for equal-mass initial black holes. For an initial infalling velocity v/c 0.462 of the colliding black 
holes, and incidence angle of collision po = 21°, we obtain a maximum Vk ~ 121 km/s located at 
n ~ 0.226. For initial equal-mass black holes (rj = 0.25) we obtain Vk ~ 107 km/s. Based on the 
integrated Bondi-Sachs momentum conservation law we discuss a possible definition of the center- 
of-mass velocity of the binary merged system and show that - in an appropriate inertial frame - 
it approaches asymptotically the net kick velocity, which is the velocity of the remnant black hole 
in this inertial frame. For larger values of v/c we obtain substantially larger values of the net kick 
velocity, e.g., for v/c ~ 0.604 a maximum Vk ~ 610 km/s is obtained. 

PACS numbers: 04.30.Db, 04.25.dg, 04.70.Bw 



I. INTRODUCTION 

The collision and merger of two black holes is presently 
considered to be an important astrophysical configura- 
tion where processes of generation and emission of grav- 
itational waves take place (cf.pj and references therein). 
The radiative transfer involved in these processes, eval- 
uated in the full nonlinear regime of General Relativity, 
shows that gravitational waves extract mass, momentum 
and angular momentum of the source, and may turn out 
to be fundamental for the astrophysics of the collapse of 
stars and the formation and population of black holes in 
galaxies 

In this vein we have recently examined numeri- 
cally the gravitational wave production and related ra- 
diative processes, as mass-energy and momentum ex- 
traction by gravitational waves, in a general collision 
of two Schwarzschild black holes 5], described in the 
realm of non-axisymmetric Robinson- Trautman (RT) 
spacetimes @. Characteristic initial data that represent 
two instantaneously boosted black holes in non head-on 
collisions were constructed for the RT dynamics. The 
initial data already present a common apparent horizon 
so that the dynamics covers the post-merger phase of 
the system. The present paper focuses on the processes 
of momentum extraction by contemplating the time be- 
havior of the momentum fluxes of the gravitational wave 
emitted and the associated total impulses imparted to 



the binary merged system. Our treatment is based on 
the Bondi-Sachs energy momentum conservation laws0- 
[Toj . The use of the Bondi-Sachs formulation gives a clear 
and simple picture of the kick processes generated in the 
recoil of the the binary merged system which are related 
to the maximum impulse due to the gravitational wave 
flux in each regime characterized by the initial data. 

The kick processes and associated recoils in the merger 
of two black holes have been investigated within several 
approaches, all of them connected to binary black hole 
inspirals. Post-Newtonian approximations (cf. (Tlj and 
references therein) estimated the kick velocity accumu- 
lated during the adiabatical inspiral of the system up to 
its innermost stable circular orbit (ISCO) plus the kick 
velocity accumulated during the plunge phase (from the 
ISCO up to the horizon); the plunge contribution dom- 
inates the recoil but the evaluation of the dynamics in 
the latter stage is the most uncertain. A recent treat- 
ment using post-Newtonian calculations and incorporat- 
ing "close-limit approximation" calculations is given in 
|12| . The first fully numerical verification of the recoil in 
nonspinning black hole binaries was reported by Baker 
et al.[2j for a mass ratio ~ 0.667 of the two black holes. 
Later Gonzalez et al. [IH undertook a more complete full 
numerical relativity examination of kicks in the merger of 
nonspinning binaries by contemplating a larger parame- 
ter domain. Recently Rezzola et al.[? ] obtained an 
important injective relation between the kick velocities 
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and the effective curvature parameter of the global ap- 
parent horizon, in head-on collisions of two black holes in 
RT dynamics using initial data derived in Refs. [l5Lll6j. 
Some features of a non head-on collision are distinct from 
those of a head-on collision [15|, [l7| or from the merger of 
inspiral nonspinning binaries, mainly for large values of 
the initial mass ratio parameter and of the incidence an- 
gle of collision of the two initial colliding black holes. One 
of them is the nonzero net momentum flux carried out 
by gravitational waves in the case of equal-mass initial 
colliding black holes. Our results will be compared with 
previous estimates of black hole kicks from full numer- 
ical relativity simulations of binary black hole inspirals 
(cf. [3] and references therein). 

The presence of kick processes and the magnitude of 
the associated recoil velocities can have important conse- 
quences for astrophysical scenarios, as the evolution and 
the population of massive black holes in galaxies or in the 
intergalactic medium[l8|, [l9j]. Observational evidence of 
black hole recoils have been reported in 0, Q and refer- 
ences therein. 

RT spacetimes@, 0] are asymptotically flat solutions 
of Einstein's vacuum equations that describe the exterior 
gravitational field of a bounded system radiating gravita- 
tional waves. In a suitable coordinate system the metric 
can be expressed as 



da 2 = [\(u,6A)- — + ^r^-)du 2 + 2dudr 
i r K 1 

r 2 K 2 {u, 9, <j))(d8 2 + sin 2 6dip 2 ). 



(1) 



where 

A(u,6>, <j>) = -i- 



1_ _ (Kg sin g /K) e __}_ («% _ 
K 2 K 2 sin 9 sin 2 9\K 4 K 3 



The Einstein vacuum equations for ([T]) result in 
„ K u 1 ((\ B sm9) g A^ \ 

In the above, subscripts u, 9 and <j) denote derivatives 
with respect to u, 9 and 0, respectively, too is the only 
dimensional parameter of the geometry, which fixes the 
mass and length scales of the spacetimc. Eq. ([3]), de- 
noted RT equation, governs the dynamics of the system 
and allows to propagate the initial data K (uq, 9, 4>), given 
in the characteristic surface u = uq, for times u > Uq. For 
sufficiently regular initial data RT spacetimes exist glob- 
ally for all positive u and converge asymptotically to the 
Schwarzschild metric as u — > oo[20|. The field equations 
present two stationary solutions which will play an im- 
portant role in our discussions: (i) the Schwarzschild so- 
lution corresponding to K — Kq — const, with A = 1/Kq 



and mass Ms, 



hw 



itiqKq; and (ii) 



K{e,<t>) = 



K 



cosh 7 + (n • x) sinh 7 ' 



(4) 



where x = (sin 9 cos <fi, sin 9 sin <j>, cos 9) is the unit vector 
along an arbitrary direction x, and n = (ni,ri2, n^) is a 



constant unit vector satisfying n\ + n^ + n\ = 1. Also Kq 
and 7 are constants. We note that ^ yields A = 1/Kq, 
resulting in its stationary character. This solution can 
be interpreted Q as a boosted black hole along the axis 
determined by the unit vector n with boost parameter 
7, or equivalently, with velocity parameter v = tanh7. 
The K(9,4>) function (Q}, which depends on three pa- 
rameters, is a K-transformation of the generalized Bondi- 
Metzner- Sachs (BMS) groupfl, and represents a gen- 
eral Lorentz boost of the BMS. 

The Bondi mass function of this solution is m(9, (f>) — 
moK 3 (9,(f>), and the total mass-energy of this gravita- 
tional configuration is given by the Bondi mass 



M B = (1/4tt) 



d<f> / m(e,(p)sme d0 



— ttioKq cosh 7 = moKo/-\/ 1 — v 2 . 



(5) 



The interpretation of ((H) as a boosted black hole is rel- 
ative to the asymptotic Lorentz frame which is the rest 
frame of the black hole when 7 = 0. This asymptotic 
inertial frame will be the one we refer to in Section 2, 
in the definition of the Bondi-Sachs energy-momentum, 
and in Section 5 when the center-of-mass motion of the 
merged system is discussed. 

In the paper we use units such that G = c = 1; c is 
however restored after the definition of the kick velocity. 
Except where explicitly stated, all the numerical results 
were done for 7 = 0.5. In our computational work we 
used mo = 10 but the results are given in terms of u/vtiq. 
We should note that we can always set too = 1 in the RT 
equation ^ by the transformation u — > u = m/to-q. 



II. THE BONDI-SACHS FOUR MOMENTUM 
FOR RT SPACETIMES, THE NEWS AND THE 
INITIAL DATA 

Since RT spacetimes describe the asymptotically fiat 
exterior gravitational field of a bounded system radiating 
gravitational waves and the initial data of its dynamics 
is prescribed on null characteristic surfaces, they can be 
included in the the realm of the Bondi-Sachs 2+2 formu- 
lation of gravitational waves in General Relativity [8H10|. 
Consequently suitable expressions for physical quantities 
to be used in the description of gravitational wave emis- 
sion processes, as the Bondi-Sachs four momentum and 
its conservation laws, must be properly derived. From 
the supplementary vacuum Einstein equations Ruu = 0, 
R ue = 0, and Ru<s> = in the 2 + 2 Bondi-Sachs 
formulation [SllSl (where (U, R, O, 4>) are Bondi-Sachs co- 
ordinates), we obtain]^, |]| 



dm(u, 9, (j>) 



+ 4c 



du 

(2) COS< 



= -^| C W 2 + c< 2 > 2 



ld_ 

2 du 



3c!, 1 ' cot ( 



2^+^ + -^ 



sin^6> 



„(i) 

-<t»t> 



(6) 



where m(u,9,<f>) is the Bondi mass function and 
ci}\u,9,(f)) and Cu\u,6,<j>) are the two news functions, 
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in RT coordinates. The extra factor K in the first term 
of the second-hand-side of Eq. ([6|) comes from the trans- 
formation of the Bondi time coordinate U to the RT co- 
ordinate u, lim r ^oo = dU/du = l/K [211 . For the RT 
spacetimes JT]), the news are expressed as 22, 23] 



2 V siri 8 



1 



c^{u,9,4>) = -T^lPg^-V^ cot l 



(7) 



where we have introduced the variable V = l/K, for 
notation convenience. We remark that c4 = = c/P 
for the boosted Schwarzschild solution (j4j, as should be 
expected. Our derivation [23| of Eqs. ([7]) extends the 
work of von der Gonna and Kramer [2lj to the case of 
non-axisymmetric RT spacetimes. and agrees with Eq. 
(45) of Cornish and Micklewright[22] modulo the use of 
the RT time derivative and of the function V adopted in 
their representation of RT metric. 

The Bondi- Sachs four momentum is defined asfioj 



d(f) / m(u,0,4>) l^sxaO dO, (8) 
Jo 



where — (1, — sin6*cos0, — sin#sin</>, — cos#), relative 
to an asymptotic Lorentz frame 1 . From Eq. (|6]) the 
Bondi- Sachs four momentum conservation law follows 



dP"(u) 
du 



(9) 



where 



where P w (u) = (Pfr(u), 0, Pfr(u)), with 



P^u) = - 



Pw(u) 



4tt 



sin 2 6»cos0 K[ ci 1)2 + c ( * )2 )d8, (12) 



cos (9 sin (9 K( c^ )2 + d 2)2 )d0. (13) 



Obviously P v (u) is conserved, this fact being a conse- 
quence of P^y{u) = for all u. 

The initial data to be used was derived in Ref. Q and 
represents two instantaneously colliding Schwarzchild 
black holes in the (x, z) plane, at u = uq, 



K(uo, 



1 



\/ cosh 7 + cos 9 sinh 7 

Q 



+ 



\J cosh 7 — (cos po cos 8 + sin po sin 6 cos 0) sinh 7 



.(14) 



For a = 0, Eq. (fT4"f corresponds to a black hole boosted 
along the positive z-axis. The parameter po of the initial 
data, denoted the incidence angle of collision, defines the 
direction of the second initial black hole with respect to 
the z-axis; a is the mass ratio of the two initial colliding 
black holes and v — tanh7 their initial infalling veloc- 
ity. This data already has a common apparent horizon 
so that the evolution covers the post merger regime up 
to the final configuration, when the gravitational wave 
emission ceases. Finally we note that, in the full Bondi- 
Sachs problem, further data - the news functions - are 
needed to determine the evolution of the system. How- 
ever for the RT dynamics the news are already specified 
once the data (fl"4"|) is given, cf. Eqs. Q. 



1 

47T 



K V- 



,(1)' 



„(2)- 



sin# dO (10) 



is the net flux of energy-momentum carried out by the 
the gravitational waves emitted. We note that the term 
between square brackets in the right-hand-side of © van- 
ishes in the integrations due to the boundary conditions 
satisfied by the news, = — and = = 
at 8 = 0, 7r. The mass-energy conservation law (Eq. © 
for fj, = 0) is the Bondi mass formula. Our main interest 
in the present paper is the analysis of the linear mo- 
mentum conservation law (Eq. ([9]) for fi = x,y,z). As 
we have discussed in Q a general collision of two black 
holes in RT dynamics is planar, namely, the dynamics is 
restricted to the plane determined by the directions of 
the two initial colliding black holes. Therefore, without 
loss of generality we fix this plane as the (x, z)-plane so 
that the momentum conservation equations relevant to 
our discussion reduce to 



dP(u) 
du 



(11) 



1 The four vector actually defines the generators (l^d/dU) of the 
BMS translations in the temporal and Cartesian x, y, z directions 
of the asymptotic Lorentz frame [To| . 



III. NUMERICAL EVOLUTION 

The initial data ([14")) is evolved numerically via the RT 
equation, which is integrated using a Galerkin spectral 
method with a projection basis space in two variables 
(the spherical harmonics) [lij . The numerical method is 
accurate and highly stable to long time runs so that we 
are able to reach the final configuration of a boosted rem- 
nant black hole when the gravitational wave emission 
ceases [5J. This long time computation is demanded in 
order to determine the behavior of basic quantities nec- 
essary for the evaluation of the processes as, for instance, 
the total impulse imparted to the merged black hole by 
the gravitational radiation. Therefore we can examine 
physical phenomena in the nonlinear regime where full 
numerical relativity simulations might present some dif- 
ficulties due, for instance, to the limited computational 
domain. In this sense the results of the evolution of the 
above data may be considered as complementary to full 
numerical relativity simulations on describing the late 
post-merger regime. 

In our numerical experiments in the present paper we 
vary a in the interval [0, 1] with fixed po = 21° and 
7 = 0.5 (corresponding to an initial infalling velocity 
v ~ 0.462). For comparison purposes, results for other 
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values of the initial parameter po are also considered. 
The truncation of the method adopted in our computa- 
tion is Np = 7. Exhaustive numerical experiments show 
that after a sufficiently long time u ~ Uf all modal coef- 
ficients of the Galerkin expansion become constant up to 
twelve significant digits, corresponding to a final time of 
computation Uf. At Uf the gravitational wave emission 
is considered to effectively cease. From the final constant 
modal coefficients we reconstruct K(uf,0,(j>) that, in all 
cases, can be approximated as 



K( Uf ,e,< 



cosh 7/ + (m / sin 9 cos (f> + nsf cos 9) sinh 7/ 



(15) 



With the final parameters (Kf,jf,riif,n3f) obtained 
from the final modal coefficients, we have in all cases 
that the rms error of Eq. (TTSj) is of the order of, or 
smaller than 10 -12 . The final configuration corresponds 
then to a boosted Schwarzschild black hole (cf. (01) with 
a final velocity Vf = taiiiryf along the direction deter- 
mined by nf — (nx/,0,n3/), and a final Bondi rest mass 
m Kj. In general 7/ < 7 and Kf > K - The angle 
Pf = arccos (^3/) - denoted the scattering angle - de- 
fines the direction of the remnant with respect to the 
z-axis. Within the numerical error of our computation 
we have (^i/) 2 + (^3/) 2 = 1, as expected. 

The values of the parameters of the remnant black hole 
are one of the basic results to be extracted from our nu- 
merical experiments that, together with the initial data, 
allow us to evaluate quantities which are characteristic of 
the radiative transfer processes involved in the gravita- 
tional wave emission. Also, by computing A(u, #,</>) for 
all u > uq, we can obtain the time behavior of important 
physical quantities, as for instance the total impulse im- 
parted to the merged binary system by the emission of 
gravitational waves. Our numerical results are displayed 
in Table 1 and constitute the basis of the analysis of the 
gravitational wave recoil, in the next Section. 



IV. GRAVITATIONAL WAVE MOMENTUM 
FLUXES AND THE KICK PROCESSES IN A 
NON HEAD-ON COLLISION 

We are now ready to examine the processes of mo- 
mentum extraction and the associated impulses imparted 
to the binary merged system by the emission of gravi- 
tational waves. Our starting point is the construction, 
via the numerically integrated function K(u, 9, 4>), of the 
curves of the net momentum fluxes carried out by gravi- 
tational waves. 

In Fig. [1] we display the curves of the net momentum 
fluxes P^(u) and P^,(u) for mass ratios a = 0.2, 0.5, 
and fixed 7 = 0.5 and incidence angle po = 21°. In both 
cases the curves present an initial regime of positive flux 
Pw(u) up to u — Uk when P^r(uk) — 0, with consequent 
increase of the linear momentum of the merged system 
along this direction, dP z (u)/du > 0. The value of Uk 
increases as a increases. 



The net momentum flux P^(u) < for all uq < u < 
Uf, so that the linear momentum of the system along 
the x direction always decreases. We then see that the 
system undergoes a dominant deceleration regime until 
it reaches the final configuration, the remnant boosted 
Schwarzschild black hole. 

The initial regime of positive flux P^(u) > (for uq < 
u < Uk) is not present in the case of larger initial data 
parameters a and/or pq. In the next Section we will 
discuss the role of this initial positive gravitational wave 
flux in engendering, when present, an inspiral branch in 
the motion of the center-of-mass of the merged system. 



A. The integrated fluxes of gravitational waves and 
the total impulse imparted to the merged system 

Integrating in time the conservation Eq. (jlip we find 
that 



P(u) - P(« ) = Mm), 



(16) 



where Iw(u) — (Iw( u )> 0> ^w( u )) ^ s the impulse im- 
parted to the binary merged system due to the momen- 
tum carried out by the gravitational waves emitted up to 
the time u. In the above 



1 f U /*2tt 

iw(w) = — / du I dd> x 

47r Ju a Jo 



sm a cos < 



K[^+^r)d9, (17) 



,(2) 



and 



1 

4tt 



du I d(f> x 
'0 

• J 1 - 



cos 9 sin 9 K\ c k V" +c (2 , ) " \d0 



(18) 



The net total impulse imparted to the system has a dom- 
inant contribution from the deceleration regimes (where 
Pyy[u) < and P^(u) < 0) and will correspond to a net 
kick on the merged system. This physically important 
behavior has already been reported in accurate numeri- 
cal relativity simulations of the merger of unequal-mass 
binary black hole inspirals0, [l2l [l3j. 



The behavior of the impulse is illustrated in Fig. [2] 
where we plot Iyy( u ) an d Iw( u ) f° r a = ^-2, Po = 21° 
and 7 = 0.5, corresponding to the time integration of 
the dashed curves in Fig. [TJ As expected the curve for 
Iyy{u) presents a local maximum at u — Uk, when the 
area measured below the curve of Pyy( u ) starts to give 
a negative contribution to the impulse, until u = u m 
when Iy^{u m ) = 0. For large u ~ Uf the curves tend 
to a constant negative value (a plateau), corresponding 
to the final configuration when the gravitational wave 
emission ceases. It is worth mentioning that Fig. [2] also 
illustrates the dominance of a deceleration regime in the 
processes. The plateau is considered to be reached when 
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TABLE I: Summary of our numerical results corresponding to an infalling velocity v/c ~ 0.4621 (7 = 0.5) and incidence angle 
of collision p = 21°. The best fit of the points (V k ,r/) to the formula (J24J) corresponds to A ~ 3.63538199, B ~ 2.46152380 
and C ~ 0.91379025, (cf. Fig. @}. 



a 


V 


Kf 


Vin/C 


v f /c 


SV/c 


Pf 


-Iw(u f )/m 


-v? 


-Iw(uf)/m 


-n 


V k (km/s) 


0.025 


0.0238 


1.045113 


0.446421 


0.446426 


1.5271 


0.37° 


1.2792 x 10~° 


0.3362 


8.1563 x 10" B 


2.1435 


2.17 


0.050 


0.0453 


1.091454 


0.430848 


0.430867 


5.5605 


0.77° 


5.6819 x 10 _b 


1.3110 


3.3361 x 10" b 


7.6973 


7.81 


0.100 


0.0826 


1.187836 


0.400195 


0.400255 


17.9815 


1.65° 


2.7635 x 10~ b 


4.9467 


1.3834 x 10~ 4 


24.7624 


25.25 


0.150 


0.1134 


1.289162 


0.370369 


0.370480 


34.0082 


2.64° 


7.4349 x 10" s 


10.4105 


3.1932 x 10" 4 


44.7116 


45.91 


0.200 


0.1388 


1.395446 


0.341541 


0.341701 


49.3652 


3.77° 


1.5575 x 10" 4 


17.1953 


5.7681 x 10~ 4 


63.6818 


65.96 


0.250 


0.1600 


1.506701 


0.313841 


0.314042 


62.6056 


5.06° 


2.8314 x 10" 4 


24.8334 


9.0763 x 10~ 4 


79.6062 


83.39 


0.300 


0.1775 


1.622933 


0.287368 


0.287600 


72.7560 


6.53° 


4.6910 x 10" 4 


32.9220 


1.3051 x 10~ J 


91.5906 


97.33 


0.400 


0.2040 


1.870356 


0.238376 


0.238634 


82.7989 


10.17° 


1.0733 x 10" a 


49.2111 


2.2569 x 10 _tl 


103.4813 


114.59 


0.500 


0.2222 


2.137755 


0.194976 


0.195219 


80.6673 


15.06° 


2.0929 x 10" a 


64.2673 


3.3174 x 10~ a 


101.8707 


120.45 


0.525 


0.2257 


2.207669 


0.185049 


0.185282 


79.1673 


16.53° 


2.4287 x 10" J 


67.71537 


3.5817 x 10~ J 


99.8653 


120.66 


0.600 


0.2344 


2.425151 


0.157534 


0.157731 


69.8608 


21.77° 


3.6700 x 10~ a 


77.1919 


4.3295 x 10~ a 


91.0621 


119.38 


0.700 


0.2422 


2.732560 


0.126622 


0.126756 


53.9660 


31.14° 


5.9634 x 10" a 


87.6814 


5.0986 x 10~ a 


74.9656 


115.36 


0.800 


0.2469 


3.059989 


0.103262 


0.103331 


35.8323 


44.09° 


9.1477 x 10~ 3 


95.7801 


5.3968 x 10" J 


56.5069 


111.21 


0.900 


0.2493 


3.407441 


0.088861 


0.088880 


1.5271 


60.85° 


1.34131 x 10~ 2 


101.7102 


4.9651 x 10~ a 


37.6501 


108.45 


1.000 


0.2500 


3.774920 


0.084214 


0.084214 


0.0016 


79.50° 


1.8965 x 10~ z 


105.7664 


3.5148 x lQ- 6 


19.6022 


107.57 



0.0006 



0.0004 



0.0002 



0.0000 



-0.0002 



-0.0004 



-0.0006 




-0.0002 



-0.0004 



-0.0006 



u/m 




FIG. 1: Linear-log plot of the net fluxes of momentum P^/(u) 
and P{jy(u) for a — 0.2, 0.5, and fixed 7 = 0.5 and incidence 
angle po = 21°. For both cases, an initial regime of positive 
momentum flux along the z axis is present, corresponding 
to an increase of the z-component of the linear momentum, 
dP z (u)/du > 0, up to u — Uk when P^(uk) = 0. The domi- 
nant contribution to the net impulse imparted to the system 
comes from the deceleration phases and will correspond to a 
net kick on the merged system. 



FIG. 2: Plot of the impulses of the gravitational waves Iw( u ) 
(dashed curve) and Iw(u) (continuous curve), for a = 0.2 and 
fixed initial incidence angle and infalling velocity, po = 21° 
and v ~ 0.462, respectively. The curve of Iyy{ u ) presents a 
local maximum at u = u k , when the deceleration regime along 
z starts. At u — u m , /fffiim) = 0. For large u ~ Uf both 
curves tend to a constant negative value, corresponding to 
the final configuration when the gravitational wave emission 
ceases. 



— Iw(u + h)\ < 10~ 10 , where h is the stepsize 
of the integration used for the evaluation of Iw(u)- At 
this stage the remnant black hole has a momentum P = 
(ni/,0,n 3/ ) Pf, with 

P f = V(P*(u/)) 2 + (P*(u/)) 2 = m Kj sinh 7/ , (19) 

whose distribution as function of a, for several po, is 
given in Q. 



From Eq. (|16[) we can derive that 

P*(u f )-P x (u ) = r w {u s ), (20) 

P*(u f )-P*(u ) = F w {u } ), (21) 

where the right-hand sides of pp]) and pi) are the 
nonzero components of the net total impulse l\y(uf) gen- 
erated by the gravitational waves emitted. 
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We remark that for larger values of a, the final im- 
pulse Iy^{uf) is larger than Iyy(uf), contrary to the cases 
a = 0.2, 0.5 shown in Fig. [21 This is shown in Figs. 
[3] where we plot the net momentum fluxes P^r (u) and 
Pw( u )' an d the associated impulses (integrated fluxes) 
Iy/{u) and 7^ (it), for the case of equal- mass initial col- 
liding black holes (a = 1) and incidence angle po = 21°. 

The purpose of Figs. [3]is three-fold. First to show that, 
for high values of the mass ratio, the initial phase of pos- 
itive gravitational wave flux along the z-axis is absent, 
with an overall deceleration of the binary merged system 
during the whole regime of gravitational wave emission. 
Second, the total final impulse along the z axis is about 
one order of magnitude smaller than the total impulse 
along the x axis; both of them are at least one order 
of magnitude larger that the corresponding impulses for 
a = 0.2. Third, in a non-head-on collision of two equal- 
mass colliding black holes, the net gravitational wave 
flux, and the associated impulses are nonzero, contrary 
to the cases of head-on collisions or inspiral binaries of 
equal-mass black holes. The implications of the third 
point will be detailed in the next subsection. 



B. The kick velocity and the modified Fitchett 
distribution for non head-on collisions 

We are now led to define the net kick velocity of 
the binary merged system as proportional to the net mo- 
mentum imparted to the system by the total impulse of 
gravitational waves up to u — Uf. We remark that our 
definition of the net kick velocity is based on the impulse 
function lyy(u) of the gravitational wave emitted evalu- 
ated at u = Uf (cf. Eqs. ([2"U|) - ([2T]) '). These definitions are 
in accordance with [13|. We obtain (restoring universal 
constants) 



V k . = 



mo if? 



lw(uf), 



(22) 



with modulus 



^I^UfY+I^UfY , (23) 



where m^K^ is the rest mass of the remnant black 
hole. We note that I w (u ) = 0. The velocity ([23]) 
is directed along an axis making the angle Of — 
arctan(/^(u/)//^(u/)) with the negative z-axis. The 
interpretation of this angle, as defining the asymptotic 
direction of the center-of mass motion in a particular 
inertial frame, will be given in Section V. From our nu- 
merical results we evaluate Vk for several values of a, 
and for fixed 7 = 0.5 and incidence angle po = 21°. The 
results are summarized in Table f where, in accordance 
with the literature, we use the symmetric mass parameter 
T] = a/(l + a) 2 . 

In Fig. j4| we plot the points (Vk,rj) from Table 1. 
The continuous curve is the least-square-fit of the points 
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FIG. 3: Linear-log plot of the net fluxes of momentum P^,(u) 
and Pw (u) (top) and the associated total impulses I^y (it) and 
Iw( u ) (bottom), for the case of equal- mass initial colliding 
black holes (a = 1) and incidence angle po = 21°. The kick 
on the merged system for this case is Vk ~ 107 km/s. 



to the analytical formula 

V = Arf{l - 4Cr7) 1/2 (l + Brf) x 10 3 km/s, (24) 

with best fit parameters A ~ 3.63538199, B ~ 
2.46152380 and C ~ 0.91379025. The maximum net an- 
tikick obtained is ~ 120.66 km/s at 77 ~ 0.2257. Eq. 
(|24l) is a modification of an empirical formula originally 
derived from post-Newtonian analysis[25j and used by a 
number of authors pdl - [l3t to describe the distribution of 
kick velocities. The additional parameter C is necessary 
to account for the nonzero net gravitational wave flux 
in non-head-on collision of two equal-mass black holes. 
Actually this modification introduced in (pH)) is the only 
one that works to produce an accurate fit, with rms error 
~ 0.0376 (and a normalized rms error ranging from 0.5% 
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FIG. 4: Plot of the points (14,77), where 14 is the net kick 
velocity due to the impulse imparted to the merged system by 
the gravitational waves, for po = 21°. The continuous curve 
is the least-square-fit of the points to the analytical formula 
flH| , with best fit parameters A ~ 3.635382, B ~ 2.461524 
and C ~ 0.913790 (cf. Table 1). The maximum of the curve 
corresponds to V k ~ 120.66 km/s for 77 ~ 0.2257 (a ~ 0.5249) 




0.05 0.1 0.15 0.2 0.25 



FIG. 5: Comparison between distributions of the net kick 
velocities 14 for a collision with incidence angle po = 21° 
(squares) and a head-on collision po = 0° (circles). The head- 
on collision is well fitted by a Fitchett 77-scaling law (C = 1) 
while the po = 21° collision is fitted by the modified Fitch- 
ett law with C ~ 0.9138 which accounts for the nonzero net 
momentum flux for equal-mass black holes. 



to 0.031%). 

For comparison purposes we plot in Fig. [5] the points 
(14,77) for the case of a head-on collision po = 0°[l7|. 
together with those of the non-head-on collision with 
po = 21° given in Fig. [4] The continuous curves are the 
best fit of the points to the analytical curve (|24p with the 
best fit parameters C — 1 for (po = 0°) (a genuine Fitch- 
ett curve) and C ~ 0.91379025 for p = 21°. As p -> 



the parameter C — > 1 monotonically, as verified numeri- 
cally. For the head-on case the distribution shows a con- 
sistent similarity with the results of Refs. [HI, [l3[ where 
kicks in black hole inspiral binaries are estimated from 
full numerical relativity simulations and from combining 
numerical relativity with perturbation theory. This con- 
sistency is mainly connected to the order of magnitude 
of the kicks, and the location in 77 of the maximum kick. 
Using a Newtonian reasoning we may conclude that the 
momentum of each body in the inspiral binary lies along 
the axis connecting the bodies and have the opposite di- 
rection, as seen from the center-of-mass frame of the sys- 
tem. On this basis we might expect that the momentum 
fluxes in this system and in a head-on collision present 
similar features. 

We should note that, since our choice of 7 = 0.5 in 
the present paper was not fixed by any physical con- 
sideration, a larger 7 in the initial data could be cho- 
sen so that the numerical values of the maxima agree 
with, or have substantially larger values than the ones 
in [13|. For instance, for 7 = 0.7 (v/c ~ 0.6044) and 
a value of a = 0.524979 (the latter corresponding to 
the location of the maximum in the case 7 = 0.5, cf. 
Fig. [2} we obtain 14 — 610.06 km/s, with associated 
total net impulses I^(u f )/m ~ -0.009753504379 and 
I^{uf)/m -0.016404424812. 

A remark is in order now, concerning the balance be- 
tween the total rest mass of the remnant and the to- 
tal net impulse of the gravitational waves in the dis- 
tribution of the net kick velocity shown in Fig. As 
can be seen from Table 1, as 77 increases from to 
0.25, both the parameter Kf and the total net impulse 
Iw{uf) = \Aw( u /) 2 + I^riuf) 2 increase. However the 
increase of the rescaled Bondi rest mass of the rem- 
nant, K% is smaller than the increase of Iw{uf) up to 
77 ~ 0.225, implying that in this range the net kick veloc- 
ity increases in accordance with (|23|) . Beyond this point 
the increase of Kj is larger than the increase of the total 
net impulse leading to a decrease in the values of 14 up 
to 77 = 0.25. 

Finally we must comment that the parameter C in 
the 77-scaling law (|2"4")1 . which assumes the value C = 1 
for head-on collisions (po = 0°) and for merging black 
hole inspirals, decreases monotonically as po increases. 
In the case p = 21° we have C ~ 0.91379025 as a 
best fit value while for the limiting case po — 90° we 
obtain C ~ 0.21165880. In fact the quantity (1 - C), 
which increases (from zero) as the incidence angle po 
increases, has a correspondence with the net kick ve- 
locity 14, or equivalently with the nonzero total gravi- 
tational wave impulse Iw(u/), for equal-mass colliding 
black holes (77 — 0.25). We mention that for the case 
pa = 90° the 77-scaling law (JUJ) with C ~ 0.21165880 
(and the other associated best fit values A ~ 1.94862882 
and B ~ 0.76602543) fits nicely the numerical points 
(77, 14) obtained. However this curve increases monoton- 
ically with 77, a behavior analogous to the momentum 
distribution of the remnant black hole for po = 90° 0]. 



8 



These results will be the subject of a future publication. 



V. THE CENTER-OF-MASS MOTION OF THE 
MERGED BINARY SYSTEM 

In our numerical work and its interpretation we were 
led to deal with several types of velocity variables, as the 
net kick velocity, the initial infalling velocity of each in- 
dividual black hole, v — tanh7, and the final velocity 
of the remnant vj = (sinp^, 0, cos pf) vt. Here we will 
further introduce the definition of an initial velocity Vi„ 
of the binary merged system when the interaction of the 
two black holes and the initial gravitational wave con- 
tent are taken into consideration. Our main interest in 
the present Section is to examine the relation among all 
these variables and try to find an appropriate definition 
of an approximate velocity of the center-of-mass of the 
merged system. The main problem in this task is con- 
nected to the fact that the system loses mass due to the 
emission of gravitational waves, which will lead us to look 
for a definition based on the momentum of the radiation 
emitted and the associated conservation laws. 

In Table 1 we have given the final velocity v/ for sev- 
eral 77, which is directed along the axis making an angle 
Pf with the positive z-axis. By using Eqs. (l2T)l) - (f2"T]) 
which are the components of the Bondi-Sachs integrated 
conservation law (|16[) evaluated at u — Uf - we can relate 
v/ to the net kick velocity through the expression 

^V t = P( U/ )-PW, (25) 

where P(uo) is the initial Bondi-Sachs momentum of the 
system which depends on a, 7 and po, and P(u/) is the 
final momentum which depends on 7^, Kf and pt. The 
final velocity associated with the final momentum P(ut) 
is given obviously by V/. 

A tentative definition of the initial velocity of the 
merged system can be calculated from the initial Bondi- 
Sachs momentum P(«o) as Vi„ = P(uo)/Mb(uo), 
where Mb(uq) is the initial Bondi mass of the sys- 
tem. P(ito) can be evaluated (with radiative correc- 
tions included) from the conservation law (fT6|) by tak- 
ing u — Uf. The modulus of Vj„ is obviously given by 
Vin = ^P x {u ) 2 + P z (u ) 2 /M B (u ). The values of v in , for 
several 77, are listed in Table 1. We note that the initial 
velocity Vi n is distinct from the previously defined ini- 
tial infalling velocity (v = tanh.7) in the following sense: 
strictly speaking the interpretation of v should actually 
be the initial infalling velocity of each hole at infinity 
before being brought into interaction, when such quan- 
tity no longer has a separately defined meaning. On the 
contrary, Vi n is the initial velocity at uq when the in- 
teraction of the two black holes and the initial gravita- 
tional wave content are already set up in the initial data. 
For illustration, for 7 = 0.5 the initial infalling velocity 
is v/c ~ 0.4621 while, for 7 = 0.5 and a — 0.3, say, 
we obtain t) ra /c ~ 0.2874 (corresponding to an effective 
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FIG. 6: Plots of the points SV = |v; n — v/| (diamonds) and 
SVcc = (Co Vin — v/) (squares) versus the parameter 77. The 
dotted line and dashed line are the least-square-fits of SV and 
SVcc to the formula (I24|l . respectively, with best fit param- 
eters (A ~ 2.20190452, B ~ 5.42498354 and C = 1) and 
(A ~ 3.08913213, B ~ 3.47815851 and C ~ 0.91126685). 
The maximum net velocity SVcc obtained is ~ 118.34 km/s 
at n ~ 0.2277. The points (black circles) fitted by the solid 
line correspond to the net kick velocity of Fig. [4] showing a 
close approximation to the net velocity SVcc- 



7i„ ~ 0.2957), Vi n /c being of the order of magnitude of 

A detailed examination of the numerical results in Ta- 
ble 1 shows that both vf and Vi n have a monotonic de- 
crease with r\ (or equivalently with a). Some features 
of the momentum extraction from the merged black hole 
by gravitational waves - connected, for instance, with the 
asymmetry of the radiation emitted (namely, with a) and 
which are present in Pf and Vk - are erased in the dis- 
tribution of Vf. This is not the case of the difference of 
velocities defined by <5V = (v/ — Vj n ), with magnitude 

SV = ^„+f 2V..-V,, (26) 

where v m • v/ = v in VfCOs(p in - pf), and p in = 
&rct&n(P x (uo)/P z (uo)). The points (6V,rj) have a dis- 
tribution which is accurately fitted by the 77-scaling law 
(f2"4")l with C = 1 (cf. Fig. [TJ]), with a maximum 
SV ~ 83.11 km/s at 77 ~ 0.2087. 

This velocity variable SV cannot however be related 
to the net kick velocity due to the following problems. 
Not only SV is much smaller than the net kick velocity 
(its maximum value being approximately two thirds of 
the maximum value of the net kick velocity) but also its 
value is approximately zero for a = 1.0 while the net kick 
velocity has a large value Vk — 107 km/s at 77 = 0.25. 
Furthermore Vi n < Vf for all 77 (cf. Table 1), which is 
not physically satisfactory since the deceleration regime 
is dominant in the post-merger phase for all initial data 
parameters. The reason for these discrepancies is that 



9 



the velocity variable <5V does not satisfy any momentum 
balance, while the net kick velocity was evaluated from 
the momentum conservation law (|16j) . 

Alternative choices that do not present such problems 
may be obtained from (| 16[) evaluated at u = Uf, which 
we express as 



v ir , 



v //Co 



C I W (Uf) 

M s (m ) ' 



(27) 



where Co = MB(uo)/(moKf cosh 7/). 

Since in the nonlinear regime mass-energy is extracted 
by the emission of gravitational radiation, we have ob- 
viously that the mass ratio Co > E what - in view of 
(I27|) - makes the variable SV inappropriate for a descrip- 
tion of the dynamics (an exception being the nonphysical 
scenario in which only momentum, and not mass, is ex- 
tracted in the emission). We are then led to consider 
appropriate net velocities constructed with v in and v/ 
as 



(5Vc = (yi n - v//Co) 



or 



SVcc = (Co v, 



v /) = 



C I W (uf) 

Mb{u ) 1 



C lwiUf) 

moKf cosh 7/ ' 



(28) 



(29) 



where the mass ratio parameter Co corrects the definition 
SV for the loss of mass in the process. 

The plots of the magnitudes SVcc and, for compar- 
ison, those of the net kick velocity Vk and of SV are 
exhibited in Fig. [B] The points of SVc are not in- 
cluded since the relative difference of its numerical val- 
ues to the ones of SVcc does not exceed 0.63% for all 
< 77 < 0.25, so the plots would not be distinguish- 
able. The dashed curve is the least-square-fit of the 
points {SVccv) (squares) to the formula (|24|) . with best 
fit parameters A ~ 3.08913213, B ~ 3.47815851 and 
C = 0.91126685. The maximum of SVcc obtained from 
the fit is ~ 118.34 km/s at 77 ~ 0.2277. This distribu- 
tion shows a close agreement with the distribution of the 
net kick velocity (black circles and the solid curve). Fi- 
nally the dotted curve is the least-square-fit of the points 
(SV, 77) (diamonds) to the law (gj) with C = 1. 

From Eq. (|25p we can then obtain an analytical rela- 
tion between the net kick velocity and SVcc and/or SVc, 
given by 



V fc = 



cosh 7/ SVcc = — Co cosh 7/ <5Vc. 



(30) 



It is worth commenting here that the close agreement of 
SVcc and/or SVc with the distribution of the net kick 
velocity, as shown in Fig. [51 is due to Co being just slightly 
larger than 1 for < 77 < 0.25. For instance, for 77 = 
0.0826 (a = 0.1) we have Co - 1.00033978 and for a = 0.5 
we have Co — 1-00317349. This is sufficient to make SVc 
and/or SVcc positive and a few kilometers per second 
smaller than the net kick velocity. We also note that 
cosh 7/ ~ 1 for the relevant domain of 77 (7 = 0.5 and 
po = 21° fixed). 

From the above results we are now in a position to sug- 
gest a possible definition for the center-of-mass velocity 



of the merged binary. Our starting point is the fact that 
the velocity of the center-of-mass of the remnant black 
hole is obviously v/ , as measured by an inertial observer 
at rest at infinity (cf. comment below Eq. [5]). Also, 
by construction, the initial center-of-mass velocity of the 
merged binary is approximately Vj n and, for reasons dis- 
cussed above, we adopt it as Co v in . Relative to this same 
initial observer at infinity, and taking into account (|27p . 
we are led to introduce 



i(ti) = Co v in — 5Vcc(u) = Co v in + 



c lw(u) 
moKf cosh 7/ 



(31) 



as a tentative definition of the center-of-mass velocity at 
any time u > uq. For u — uj and u = uq we have by 
construction that v cm (ity) = v* and v(ii ) = Co v m ~ 
Vj n . According to this definition, for an inertial frame 
£ with velocity Co v m relative to a rest inertial frame at 
infinity, the motion of the center-of-mass is given by 



i(u) ~ -SVcc(u) 



c Iw(u) 
moKj cosh 7/ 



(32) 



The final velocity of the remnant (in this inertial frame 
C) is then 



Vcm(ll/) = v/ ~ - SV C c(uf) ~ V fc 



(33) 



namely, the net kick velocity, as expected 2 . 

We are now in a position to describe the motion of 
the center-of-mass of the merged binary system in the 
inertial frame C. This inertial frame is actually the rest 
frame of the center of mass at u — uq, and the trajectory 
of the center-of-mass is an integral curve of the velocity 
field 



c lw{u) 
moKj cosh 7/ 



(34) 



An integral curve of the velocity field (|3"4")l with ini- 
tial conditions x cra = = z cm is illustrated in Fig. 
[7] for initial data parameters a = 0.2 and po = 21°. 
This velocity field is proportional to the field c lw{ u ) = 
c (Iw(u), 0, I^y(u)) so that the analysis of its integral 
curve can be discussed from the behavior of the curves 
in Fig. [21 The positive impulse along the z-axis from 
u = uq to u = u m (resulting from the initial phase of 
positive gravitational wave flux along z) is responsible 
for the inspiral branch in the positive z semi-plane. For 
u — > Uf the curve approaches the asymptote which makes 
an angle 9/ = arctan(/^ / (7t/)//^ / (uf)) ~ 15.11° with 
the negative z-axis of the inertial frame C. This asymp- 
tote is the direction of the net kick velocity, which is the 
velocity of the remnant in the inertial frame C. 



2 We remark that the approximate equality signs in (132 P — (T34l) come 
from the approximations made in the relativistic addition of ve- 
locities. 
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FIG. 7: Plot of the trajectory of the center-of-mass of the 
binary merged system (a = 0.2 and po = 21°) in the inertial 
frame which is the rest frame of the center of mass at u = 
uo. The trajectory is an integral curve of the velocity field 
Vcm(w), Eq. p4[) . with initial conditions x C m = = z C m. The 
initial phase of positive impulse along the z-axis is responsible 
for the initial inspiral branch of the trajectory in the z > 
semi-plane. The trajectory has an asymptote as u — > Uf, 
which makes an angle O/ = arctan(7~^(ii/)//^(ii/)) with 
the negative z-axis of the inertial frame C (cf. text). In the 
present case G/ ~ 15.11°. 



VI. FINAL DISCUSSIONS AND CONCLUSIONS 

The present paper extends and completes [5] by a de- 
tailed examination of the gravitational wave recoil and 
the associated momentum extraction in the post-merger 
phase of two black holes in non-head-on collision, in the 
realm of RT spacetimes. The characteristic initial data 
for the system contain three independent parameters, the 
mass ratio a of the two initial colliding black holes, the 
initial infalling velocity parameter 7 and the incidence 
angle po. Our analysis is based on the Bondi-Sachs mo- 
mentum conservation laws which give a clear picture of 
the processes that generate kicks in the binary merged 
system. The impulse imparted to the merged system by 
the net momentum fluxes carried out by gravitational 
waves are evaluated as a function of time, for a large do- 
main of the initial data parameters. By numerically inte- 
grating the non-axisymmetric RT equation (via a numer- 
ical code based on the Galerkin spectral method) we eval- 
uate the total net impulse at the final time Uf, when the 
gravitational wave emission ceases. Typically the total 
net impulse is negative what, in terms of the Bondi-Sachs 
momentum conservation law, corresponds to a dominant 
deceleration regime of the system due to the emission of 
gravitational waves. However, for relatively small values 
of the parameters a and po an initial positive impulse in 
the z direction may be present, which will be responsible 
for an initial inspiral branch in the motion of the center- 
of-mass of the system. A novel feature of non-head-on 



collisions (po 7^ 0°) is the nonzero net gravitational wave 
fluxes for equal-mass colliding black holes, contrary to 
the cases of head-on collisions and black hole inspiral bi- 
naries. 

Based on the total net impulse we evaluated the net 
kick velocity[l3| imparted to the system by the gravita- 
tional waves emitted. The distributions of the net kick 
velocity V& is obtained as a function of the symmetric 
mass parameter 77 and are well fitted by a post-Newtonian 
motivated formula [l2l [25j | with an appropriate modifica- 
tion to account for the nonzero gravitational wave mo- 
mentum flux for equal-mass colliding black holes in non- 
head-on collisions. The modified Fitchett 77-scaling law 
adopted produces an accurate fit of the points (^,77), 
with a rms error ~ 0.038. The maximum net kick ob- 
tained was Vk ~ 121 km/s located at 77 ~ 0.2257, while 
for initial equal-mass black holes (77 = 0.25) we obtained 
Vk ~ 107 km/s, for an incidence angle po — 21°. Con- 
cerning the order of magnitude of the velocities and the 
location of the maximum i n 77 , these values are consistent 
with results of Refs. [HI, [l3| which treated the case of 
gravitational wave recoil in the merger of black hole bi- 
nary inspirals. The larger values of the kick velocities ob- 
tained in these references is possibly due to the fact that 
we have restricted our numerical integration to a fixed 
initial data parameter 7 = 0.5 (or an initial infalling ve- 
locity v/c ~ 0.462). However we give numerical evidence 
that for larger 7's the maximum kick velocity increases 
substantially. In fact, for 7 = 0.7 and a = 0.524979 (the 
latter corresponding to the location of the maximum in 
the case 7 = 0.5, cf. Fig. |4|) we obtain Vk — 610.06 km/s. 

The parameter C in the modified Fitchett formula 
which assumes the value C — 1 for head-on colli- 
sions and binary black hole inspirals, decreases monoton- 
ically as the incidence angle po increases. Specifically the 
quantity (1 — C), which increases (from zero) as the inci- 
dence angle po increases, has a correspondence with the 
net kick velocity Vk of equal-mass colliding black holes 
(rj = 0.25). We mention that, for the limit case po = 90°, 
the points (V k ,rj) are well fitted by (JMJ) with C ~ 0.2117 
but the fit curve is monotonic in 77, a behavior analo- 
gous to the momentum Pf of the black hole remnant. A 
complete and detailed analysis of the kick velocity distri- 
bution for a large domain of the initial data parameter 
space, with special focus on its dependence on the inci- 
dence angle po, is currently under investigation. . 

Also based on the Bondi-Sachs momentum conserva- 
tion law we have introduced a tentative definition of the 
center-of-mass velocity of the binary merged system. In a 
inertial frame that closely approximates the zero initial- 
Bondi-momcntum frame, the center-of-mass velocity (i) 
is proportional to the gravitational wave impulse vector 
Iyy(u) (or the integrated momentum flux vector Pw(u))', 
and (ii) approaches the net kick velocity as u —¥ ut . This 
inertial frame is actually the rest frame of the center- 
of-mass of the merged system at u — uq and the tra- 
jectory of the center-of-mass is an integral curve of the 
velocity field v cm = dw / (moKj cosh jf). This integral 
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curve exhibits an initial inspiral branch in the positive z 
semi-plane whenever an initial phase with P^(u) > is 
present. 

We must do a final comment on the connection between 
the simulations of the present paper (post-common- 
horizon, head-on and non-head-on collisions in RT dy- 
namics) and generic black hole merger kicks. RT evolu- 
tions do not in general model the same physical situation 
as numerical relativity (NR) evolutions of merging binary 
inspirals, and do not in general cover the same compu- 
tational domain. In fact, although the physical process 
that generates recoils in RT evolutions is the same as the 
one that generates recoils in NR evolutions, namely, the 
balance between the variation of the momentum of the 
system and the momentum flux of the gravitational wave 
emitted, the generation of the momentum flux of the 
wave may be completely distinct in each case, depend- 
ing on the specification of the required initial data for 
the dynamics. However the generated recoils in both ap- 
proaches, specifically head-on collision of two black holes 
(RT) and merging binary inspirals (NR), have impres- 
sive similar features. This suggests that - once a common 
horizon is formed - the generation of the momentum flux 
of the wave (which would determine the pattern of the re- 
coils) might present in this domain a generic behavior and 



that the characteristic approach can give complementary 
and reliable results on the momentum extraction of the 
system. Also, as discussed in the present paper, char- 
acteristic initial data for non- head-on collisions [j| satis- 
fies a slightly modified Fitchett 77-scaling law (l24l) . the 
modification being necessary to account for the nonzero 
net momentum flux for equal-mass initial colliding black 
holes; the distribution approaches the original Fitchett 
law as the incidence angle of collision approaches zero 
(the case of a head-on collision). Recently the resulting 
recoil velocity associated with a family of parametrized 
initial data (corresponding to a combination of the head- 
on collision characteristic initial data) was examined [26j 
and was shown not to satisfy Fitchett 's 77-scaling law, al- 
though leading to a zero final recoil for equal-mass black 
holes. 
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